Nonlocal correlations in normal- metal superconducting systems 
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We examine nonlocal effects between normal- metal gold probes connected by superconducting 
aluminum. For highly transparent Au/Al interfaces, we find nonlocal voltages that obey a spatial 
and temperature evolution distinct from the nonequilibrium charge imbalance signals usually found 
in such systems. These voltages are consistent with the predicted effects of crossed Andreev reflec- 
tion and elastic cotunneling, effects that involve coherent correlations between spatially separated 
electrons. 

PACS numbers: 74.45. +c, 03.67.Mn, 74.78. Na 



Nonlocal, coherent effects are hallmarks of quantum 
mechanics. Recently, interest has grown in nonlocal pro- 
cesses that coherently couple electrons in spatially sep- 
arated normal metals (N) linked by a superconductor 
(S).[100B| The correlations induced between the spa- 
tially separated electrons are predicted to give rise to 
nonlocal voltages in response to a drive current below 
the superconducting transition. Here we present exper- 
imental evidence for such nonlocal voltages in normal- 
superconducting systems with highly transparent inter- 
faces, and demonstrate that they exhibit the predicted 
spatial evolution. 

Figure GJa) and (b) illustrate two processes that are 
predicted to coherently couple electrons in normal met- 
als connected to a common superconductor. In crossed 
Andreev reflection (CAR), electrons from one normal 
metal are coupled to electrons in the other normal metal 
through a nonlocal analog of Andreev reflection, so long 
as the width of the superconductor between the two nor- 
mal metals is comparable to or less than the supercon- 
ducting coherence length, £5. In elastic co-tunneling 
(EC), electrons with predominantly sub-gap energies tun- 
nel directly from one normal metal to another. In 
both cases, the electrons remain in a coherent entan- 
gled state, coupled through their interaction with the 
superconductor. [2j For CAR, a current driven from one 
normal metal into the superconductor induces a current 
from the second normal metal into the superconductor 
due to the coupling of pairs of electrons in the two N 
leads; for EC, the direction of the induced current is op- 
posite. In both CAR and EC, the width of the supercon- 
ductor between lead 1 and 2 must be comparable to the 
superconducting coherence length, 

Experimentally, it is more convenient to measure the 
voltage on the second lead with respect to the supercon- 
ductor: for CAR, the drive current will induce a negative 
voltage on the second N lead, while for EC, the voltage 
is positive. The voltage is nonlocal in the sense that 
no net current flows through the second NS interface or 
the second N lead in the measurement. Since CAR and 
EC depend on the spin orientation of the electrons, both 
processes are affected if there is a finite spin polariza- 



tion in the normal leads, such as would be the case if 
the normal leads were ferromagnetic (F). For example, 
in the extreme limit of half-metallic ferromagnets where 
only one spin orientation is present in each ferromagnetic 
lead, EC is favored and CAR is completely suppressed if 
the magnetizations of the two ferromagnets are paral- 
lel, while CAR is favored and EC completely suppressed 
if the magnetizations are antiparallel. For extended NS 
contacts the nonlocal voltages arising from EC and CAR 
are predicted to decay exponentially with the width of 
the superconductor, L, over the length scale of £s-5\ 

A third process may complicate the observation of 
CAR and EC in NS and FS structures. This is the phe- 
nomenon of charge imbalance (CI). 0,0,0] ^ ne P n Y s i cs °f 
this process is shown schematically in Fig. HJc). A cur- 
rent, /, of quasiparticles injected from the normal metal 
gives rise to a nonequilibrium distribution of quasiparti- 
cles and Cooper pairs in the superconductor near the 
NS interface, which in turn leads to different electro- 
chemical potentials \i qv and p cp near the NS interface 
for the quasiparticles and Cooper pairs. p qp and 

jJLcp can 

be probed by normal metal and superconducting probes 
respectively. p cp relaxes to the electrochemical potential 
in the bulk of the superconductor /io over the length scale 
£5, which is typically of the order of 100 nm for conven- 
tional superconductors. p qp , on the other hand, relaxes 
to /io over the charge imbalance length Aq*, which can 
be much longer. 0, Oil If the potential at the interface 
is taken to be zero, u qp in the superconductor has the 
spatial dependence jl l| 



l^qp(x) = eA Q * p s I t&nh(x / Aq* 



(i) 



where x is the distance from the interface and ps is the 
resistance per unit length of the superconducting wire in 
its normal state. 

From the equation above, the asymptotic value of p qp 
deep in the superconductor is eAg*psI = /io- A su- 
perconducting lead used to measure the voltage in the 
superconductor near the NS interface will probe p cp ; be- 
yond a distance from the interface, this is essentially 
/io- On the other hand, a normal voltage lead placed on 
the superconductor will measure p qp . In this case, there 
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FIG. 1: (color online) Schematic representation of electronic 
processes that contribute to nonlocal voltages in normal- 
metal/superconductor structures, (a) In crossed Andreev re- 
flection (CAR), an electron from one normal metal lead (1) 
combines with an electron of opposite spin orientation from 
the second normal lead (2) to form a Cooper pair in the su- 
perconductor, (b) In elastic cotunneling (EC), a electron 
tunnels from normal lead 1 to normal lead 2 through the 
effective tunnel barrier represented by the superconducting 
gap. (c) In charge imbalance (CI), a quasiparticle current 
injected from the normal metal into the superconductor re- 
sults in nonequilibrium electrochemical potentials \i qv and \i cv 
(for quasiparticles and Cooper pairs respectively) that relax 
to the electrochemical potential in the bulk of the supercon- 
ductor on different length scales. The schematic depicts local 
CI; a nonlocal CI contribution also exists in regions of the su- 
perconductor near the NS interface not in the current path. 



will be a difference fi qp (x) — /io between the voltage mea- 
sured in the superconductor by the normal lead and any 
superconducting lead located more than a distance £s 
from the interface. The resulting enhancement in the re- 
sistance will decay with the distance x from the interface 
as AR(x) = p s A Q *(l - tanh(x/AQ*)). 

Aq* is expected to diverge at temperatures near the 
transition temperature T c of the superconductor or at 
quasiparticle energies near the superconducting gap.Q 
The CI contribution vanishes in the limit of low temper- 
atures and low energies, as shown by a number of exper- 
iments over the last thirty years. 0, E3, EH E3, E3, E3, E3| 
In these experiments, however, the induced charge im- 
balance was local, in the sense that the nonequilibrium 
distributions were probed along the current path, as de- 
picted in Fig. [He). Intuitively, one would also expect a 
nonlocal CI contribution in parts of the superconductor 
near the NS interface but outside the current path, al- 
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FIG. 2: (a) Temperature dependent four-terminal resistance 
of a NS junction shown in the inset. For the measurement, 
a small ac current is sent from one Au lead to one Al lead 
and the resulting four-terminal resistance is measured with a 
bridge using a lock-in amplifier with the two remaining leads. 
The size bar in the inset is 300 nm. (b) Differential resistance 
dV/dl as a function of the dc current I through the junction 
measured at 300 mK. / is injected in parallel to the ac current 
for this measurement. 

though the effective Aq* might be different. A detailed 
theoretical description of nonlocal CI is yet to be worked 
out. 

Two groups have recently reported observations of 
CAR and EC in patterned metallic heterostructures. 
Beckmann et a/.[l^ measured the nonlocal resistance of 
a device consisting of multiple ferromagnetic Fe wires 
in contact with an Al wire, where the magnetization of 
the Fe wires could be aligned parallel or antiparallel by 
means of an external magnetic field. They observed a 
very small difference AR in the nonlocal resistance be- 
tween the parallel and antiparallel configurations below 
T c of the Al wire, consistent with that expected for CAR 
and EC. This signal was superposed on a CI contribu- 
tion that was almost three orders of magnitude larger. 
Russo et fli.|l7j| measured the nonlocal resistance of NS 
structures with tunneling interfaces as a function of the 
voltage bias. They observed a positive nonlocal voltage 
at zero bias that turned into a negative nonlocal volt- 
age with increasing voltage bias, which they ascribed to 
a transition from EC to CAR. This is an unexpected re- 
sult in the tunneling limit, where theory predicts that EC 
and CAR are equal and opposite in sign at all energies, 
so that no nonlocal voltage should be observed. Q 

We report here measurements on NS devices with 
highly transparent interfaces, which we have found show 
large nonlocal voltages below the superconducting tran- 
sition as a function of both temperature and applied cur- 
rent bias. Two specific types of nonlocal signals are ob- 
served: the first type is more prominent at very low tem- 
peratures and current bias and decays on a short length 
scale comparable to £5. The second type of nonlocal volt- 
age is largest near T c and near the critical current and 
decays on a longer length scale. Based on their depen- 
dence on temperature, current bias and distance, we can 
readily distinguish the two types of nonlocal signals, and 
identify the first type as arising from EC and CAR, and 
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FIG. 3: (color online) (a) Scanning electron micrograph of a 
device with six nonlocal leads, each separated by 210 nm. For 
the measurements, a current I is sent through the lower Au 
lead into the superconductor as shown, and the resulting four- 
terminal nonlocal resistance is measured between the Au leads 
marked 1 through 6 and the superconducting lead labeled V—. 
A small ac current is used for the temperature dependence 
measurements, and a dc current / is injected in parallel to 
the ac current for the differential resistance measurements. 
The size bar is 1 /mi. (b) Temperature dependence of the 
nonlocal resistance of the sample in (a) measured using the 
first four leads. The numbers labeling the curves correspond 
to the lead numbers in (a). 

the second type as arising from CI. The EC/CAR signal 
in our devices is large and positive, as predicted by the- 
ory for the case of transparent contacts, where the EC 
contribution is expected to be larger than CAR.[l8| 

The characteristic difference between the EC/CAR 
contributions and the CI contribution can be seen in the 
temperature-dependent resistance of the simple NS cross 
device shown in Fig. Et a )- Above T c , the normal state 
resistance Rn of the interface is 70 mQ. This should be 
compared to a rough estimate of the Sharvin resistance 
of the NS junction based on its dimensions, which is of 
the order of 10 mSl, indicating that the NS interface is 
highly transparent. As the junction is cooled through T c , 
a large resistance peak is observed. This resistance peak, 
due primarily to the CI contribution, decreases as the 
junction is cooled further, but then saturates at around 
1 much larger than Rn- As noted earlier, the CI con- 
tribution is expected to vanish as T — > 0, hence the low- 
temperature resistance arises from EC and CAR. Similar 
behavior is observed in the low-temperature differential 
resistance dV/dl of the junction as a function of the dc 
current /, shown in Fig. 2(b). The zero-bias differen- 
tial resistance is approximately 1 ^, reflecting the low- 
temperature limit of the resistance seen in Fig. Et a )- As 
/ is increased, a peak in dV/dl develops, corresponding 
to the CI contribution of Fig. Et a )- As / is increased fur- 
ther, however, a large negative differential resistance is 
observed before the junction returns to its normal- state 
resistance at I ~ 6 jaA. The differential resistance peak 
and dip are primarily associated with CI, although we 
cannot preclude the possibility that contributions from 
EC/CAR may be present as well. 
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FIG. 4: (color online) (a) Nonlocal differential resistance of 
the device of Fig. ffla) at 20 mK. The numbers refer to the 
Au probes in Fig. Eta) with which the signal was measured. 
Inset: expanded view of the low current bias regime, showing 
that the zero-bias differential resistance decays rapidly as one 
moves away from the NS interface through which the current 
flows, (b) Zero-bias (1=0) resistance (circles) and finite-bias 
resistance peak (triangles) from (a) as a function of the dis- 
tance, L, of the nonlocal probe from the current-carrying NS 
interface. 

In order to distinguish between EC/CAR effects and 
CI, we fabricated NS devices that enable us to measure 
the truly nonlocal differential resistance as a function of 
the distance from the NS interface. Fig. Et a ) shows 
a scanning electron microscope (SEM) image of one of 
these devices. These samples were fabricated using mul- 
tilevel e-beam lithography on oxidized Si substrates. The 
50 nm thick Au layer was patterned and deposited first. 
After further patterning for the Al layer, an in situ Ar 
plasma etch was performed prior to the deposition of the 
50 nm Al layer to ensure transparent interfaces. Figure 
Efb) shows the nonlocal resistance as a function of tem- 
perature measured using the four normal metal probes 
closest to the NS interface through which the current is 
sent. The curves are similar to the data for the NS junc- 
tion shown in Fig. Et a )- In this case, however, the nonlo- 
cal resistance vanishes above T c ~ 0.6 K, as expected in 
the normal state. (We note that T c for the superconduct- 
ing wire in this device is suppressed from its typical value 
of ~1.2 K for Al due to the inverse proximity effect aris- 
ing from the large number of normal probes in contact 
with the superconductor.) Below T c , a peak develops; as 
the temperature is lowered further, the resistance satu- 
rates at a finite value. Measurements of the differential 
resistance of the interface itself give an even larger low- 
temperature resistance similar to that seen in Fig. Eta), 
indicating that the origin of the large zero-bias offsets 
seen in Figs. Et a ) and|3tb) is the same. Both the mag- 
nitude of the resistance peak and the low-temperature 
nonlocal resistance decrease as one moves further away 
from the NS interface; however, the length scales over 
which they decay are different, with the length scale for 
the low-temperature resistance being shorter as can be 
seen from Fig. Gt D ). 

Similar behavior is observed in the nonlocal differen- 
tial resistance dV/dl as a function of the dc current / 
sent through the interface. Figure EJa) shows the result- 
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ing nonlocal differential resistance at a temperature of 20 
mK for all six voltage leads shown in Fig. |3Ja). Due to 
a small asymmetry in dV/dl vs. I that likely arises from 
thermopower effects, these curves show only the symmet- 
ric component of dV/dl. The curves are similar in form 
to the curve for the NS cross shown in Fig. Eljb), but 
show the decrease of the nonlocal differential resistance 
with distance from the NS interface. While the finite- 
bias peak is observable even to the sixth lead, the inset 
to Fig. Efa) shows that the zero-bias offset essentially 
vanishes by the third lead. The different dependence of 
these two nonlocal contributions on the distance from the 
interface, L, can be seen more explicitly by plotting the 
magnitude of the zero bias differential resistance and the 
finite-bias peak resistance as a function of L. These data 
are shown in Fig. E^b). The zero-bias resistance clearly 
decays on a much shorter length scale than the finite-bias 
peak resistance. The theoretical decay of EC/CAR de- 
pends upon the geometry of the sample; for the samples 
measured here the NS interface is comparable in its di- 
mensions to £s, approaching the extended contact limit 
where an exponential decay is predicted, p] The solid line 
shows the best fit of the zero-bias resistance to an expo- 
nential decay, e~ L ^ s , with a value of £s= 315 nm. A 
slightly better fit, yielding ^5= 217 nm may be found by 
adding a small, constant zero-bias resistance to the ex- 
ponential decay. Typical values for £s for the diffusive 
Al films with T c ~ 1.2 K deposited in our laboratory are 
140-310 nm.[ll[ The fitted values of £s are in very good 
agreement with these values given that T c =0.6 K due to 
the inverse proximity effect. 

The dotted line in Fig. E{a) shows a fit of the peak 
resistance to (1 — tanh(x/Ag*)), the form expected if 
the Au leads probe the quasiparticle CI potential ja qp (x), 
with a value of Aq*= 1085 nm. We note that an even 
better fit is obtained with a linear decay, (Aq* — x), with 
Aq*= 1442 nm. We have no explanation for this ob- 
servation at present, but we note that a linear decay is 
expected if the Au leads probe the normal-state poten- 
tial iin(%) rather than the quasiparticle potential fi qp (x). 
The values of Aq* obtained from the fits are shorter than 
the value of 10-20 fim typically quoted for Al from experi- 
ments on charge imbalance phenomena performed almost 
three decades agQ.|9l ll0l.ll4j More recent experiments on 
charge imbalance are consistent with An* ~ \ ^m, in 
agreement with the values we obtain. [ill [13 

We can therefore distinguish the EC/CAR contribu- 
tion to the nonlocal resistance from the CI contribu- 
tion based its temperature, current and scaling with 
the distance from the current-carrying NS interface. 
The EC/CAR contribution at low temperature and zero 
current-bias is positive, indicating that the magnitude of 
the EC contribution is larger than that of the CAR con- 



tribution, as predicted for transparent NS interfaces. |l8j 
We cannot preclude the fact that there may also be con- 
tributions due to EC and CAR at finite bias. Indeed, 
there are experimental indications that there are contri- 
butions other than those due to CI in the current regime 
where the resistance peak and dip are observed in the 
differential resistance. If we scale the the dV/dl curves 
so that the resistance peaks for all curves lie on top of 
one another, the resistance dips will not similarly scale, 
as would be expected if both resistance peaks and dips 
were due to CI. Thus, it is quite likely that there may 
be contributions due to EC and CAR at finite current 
bias. These contributions are not precluded by the the- 
ory, although the detailed dependence on temperature 
and current bias in the regime of transparent interfaces 
is still to be worked out. 
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